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The gamma-ray excess observed by the Fermi-LAT in the Galactic Center can 
be interpreted by the dark matter annihilation to bb via a light pseudoscalar in the 
NMSSM. It is interesting to note that the corresponding singlet scalar is useful to 
achieve a strongly first order phase transition required by the electroweak baryo- 
genesis. In this paper, we investigate the possibility that the NMSSM model can 
simultaneously accommodate these two issues. The phase transition strength can 
be characterized by the vacua energy gap at zero temperature and be sufficiently 
enhanced by the tree-level effect in the NMSSM. We find that the annihilation of 
Singlino/Higgsino DM particles occurring close to the light pseudoscalar resonance 
is favored by the galactic center excess and the observed DM relic density, and a 
resulting small «/A and a negative A, can also lead to a successful strongly first 


order electroweak phase transition. 


PACS numbers: 95.35.+d, 12.60.Jv, 11.10.Wx 


I. INTRODUCTION 


The weakly interacting massive particles (WIMPs) are the most popular and attractive 
dark matter (DM) candidates. WIMPs with masses of O(1) — O(10*) GeV naturally provide 
the correct DM relic density via the freeze-out mechanism in the early Universe. The gamma 


rays produced by present DM annihilations are very good probe to reveal the microscopic 


nature and the distribution of DM due to the simple propagation and low energy loss. Such 
DM induced gamma-ray signatures from astrophysical sources with high DM densities, such 
as dwarf galaxies, galaxy clusters and the Galactic Center, have been extensively studied. 

Recently, an extended gamma-ray excess of a few GeV in the GC was discovered in the 
Fermi Large Area Telescope data [IHI]. Although the galactic center excess (GCE) is 
highly statistically significant, the diffuse gamma-ray background model in the GC would 
introduce some unclear systematic uncertainties. After considering such uncertainties, some 
studies still confirmed the existence of the GCE [I]. The origin of the GCE may be 
astrophysical sources, such as a population of millisecond pulsars (MSP) [13]. However, 
the MSP scenario is not easily explained the energy spectrum and spatial distribution of the 
GCE [14] [15]. Additionally, whether emissions from all MSPs can account for the total flux 
of the GCE is still under debate [13]. 

A very attractive interpretation for the GCE is the DM annihilation. The spatial distri- 
bution of the GCE is compatible with the square of Navarro-Frenk-White DM distribution 
with a slop ~ y = 1.2. The energy spectrum of the GCE is well fit by a ~ 30 — 40 GeV DM 
particle annihilating into bb with a cross section of (ov) œ 2 x 107-76 em3/s. A ~ 7—10 DM 
particle annihilating into T7 with a cross section of (ov) ~ 5 x 10-2” cm/s is also possible 
to fit the GCE, but the quality of this fit is lower than the fit for bb channel [8]. These 
models can be directly tested by the Fermi-LAT observations of dwarf spheroidal galaxies. 
The latest Fermi-LAT results have improved the current limits by a factor of O(1) and set 
constraints on the DM models accounting for the GCE [16]. However, considering the un- 
certainties from the DM density profiles in the GC and dwarf galaxies, the tensions between 
different Fermi-LAT observations can be relaxed [17]. 

In literatures, simplified models containing new mediators connecting DM particles to bb 
have been proposed to explain the GCE [18H20]. Since the masses of mediators required 
by the GCE are within the region of O(107) — O(10*) GeV, they can be directly on-shell 
produced at colliders and then decay into DM particles. At direct detection experiments, 
even there is no directly interaction between the DM particle and light quarks at tree level, 


DM-b quark interactions may also induce DM-nucleon scatterings at loop level. Therefore, 


the results of collider and direct detection experiments would set constraints on simplified 
models accounting for the GCE [I8H20]. 

In principle, UV complete theories can be mapped into simplified models [21]. Super- 
symmetry (SUSY) is an very attractive and theoretical motivated new physics model; it can 
provide suitable DM candidates and mediators to explain the GCE. In the SUSY model, 
DM particles can annihilate into bb via t-channel exchange of light sbottoms. The limits 
on directly sbottom pair production from LEP and LHC results are very stringent. A more 
promising s-channel annihilation channel is through a pseudoscalar yy — a — bb. Compar- 
ing with the annihilation mediated by a scalar, this process is not suppressed by the small 
DM velocity in the Galaxy. Since the pseudoscalar mediator accounting for the GCE is 
required to be lighter than ~ 100 GeV, the LHC Higgs searches have set strong limits on 
the pseudoscalar in the context of the minimal supersymmetric standard model (MSSM). 
An economical SUSY extension model evading all collider and direct detection limits is the 
next-to-minimal supersymmetric standard model (NMSSM) [22], which provides a correct 
SM-like Higgs mass ~125 GeV and solves the u problem in the minimal supersymmetric 
standard model. In the NMSSM, the pseudoscalar mediator is the lightest CP-odd Higgs. 
As pointed in Ref. 21], the GCE and DM relic density can be simultaneously interpreted 
by the Singlino/Higgsino DM resonant annihilation or the Bino/Higgsino DM off-resonant 
annihilation via the exchange of light pseudoscalar (see also Refs. [23}26]). 

Another fundamental problem in the Cosmology is the origin of the baryon asymmetry. 
An attractive scenario is the electroweak baryogenesis (EWBG), where a strongly first order 
electroweak phase transition (SFOEWPT) is required to avoid the washout of the generated 
baryon asymmetry of Universe (BAU)! . The EWBG connects the Cosmology to the Higgs 
physics and can be tested at high energy colliders. In the MSSM, a light stop is required 
to successfully trigger the SFOEWPT [2831]. After the Higgs discovery, a big challenge 
is how to simultaneously explain the correct Higgs mass ~ 125 GeV while evading strin- 


gent constraints from Higgs measurements. In fact, recent results of LHC Higgs and stop 


1 Tn fact, SFOEWPT is not enough for the successful EWBG and the strength of the SFOEWPT may affect 
the CP violation source which drives the EWBG [27]. 


searches have ruled out the mass window of the light stop in the MSSM[32H35]. If the MSSM 
Higgs sector (consisting of two Higgs doublets) is enlarged to include an additional singlet 
superfield, we have an alternative supersymmetric framework, namely the next-to-minimal 
supersymmetric Model (NMSSM), which could provide the successful SFOEWPT, as been 
appreciated for a long while [86H43]. As pointed in Ref. [42], the SFOFWPT would occur in 
R, ~ —1 and positive R, < O(10) in NMSSM, where Rẹ, = 4kv,/A, is a critical parameter. 
In this case, the lightest CP-odd Higgs is very light. It is interesting to find the connec- 
tion between the DM phenomenology and the EWBG. In this paper, we have studied this 
problem in detail and found the suitable parameter space simultaneously accommodating 
the GCE, DM relic density and SFOEWPT in the NMSSM. 

This paper is organized as follows. In Sec. II, we reviewed and resummarized the key 
of the strong first order electroweak phase transition in the NMSSM, and we make the 
energy gap analysis which gives the clue of EWPT. In Sec. III, the interpretations of GCE 
in NMSSM is reviewed, the resonance effect in the situation of K/A < 1 is analyzed. In 
Sec. IV, the numerical analysis of the EWPT, GCE and DM relic density are carried out, the 
benchmark scenario which could explain GCE and give rise to the SFOEWPT and correct 
DM relic density is presented. In Sec.V, we summarized the work and give our discussions 


and conclusion. 


II. STRONGLY FIRST ORDER ELECTROWEAK PHASE TRANSITION IN 
THE NMSSM 


The NMSSM model can solve the u problem in the MSSM and provide a 125 GeV Higgs 
boson without large loop corrections. Throughout this paper we restrict ourselves to the Z3 


NMSSM, where an extended superpotential is given by 
1 
W = ASH,Ha + Bho, (1) 


After the singlet filed S getting a vacuum expectation value (VEV) vs, an effective u term 


can be naturally generated 


U = pape = AUVs, (2) 


which is expected to be of O(107) GeV. The soft breaking terms in the Higgs sector are given 
by 
1 
Veo: = AACA Ha + gh Ans” + h.c. (3) 


The angle £ is defined as 

v, =vsin Pp, va = v cosp, (4) 
where v, and vq are VEVs of H, and H4 respectively, and v = „/v2 +v? = 174 GeV. 
Compared with the MSSM, the tree-level Higgs mass in the NMSSM is enhanced by a new 


term \”|H°H9|?. Therefore, a 125 GeV Higgs can be easily achieved in the NMSSM. The 
mass square of the SM-like Higgs is given by 


mi, = mz cos” 28 + Av? sin? 28 + Smog + OM ins (5) 


where ÔM op and MZ ix denote the loop effects and the mixing effects on the Higgs mass. 
The dynamics of the EWPT are governed by the finite temperature effective Higgs po- 
tential, which reads 


Vert = Vitee + Vow + Vor + Vr + Vaaisy- (6) 


In Eq. g Vīree denotes the tree-level Higgs potential in the NMSSM. Vew is the well-known 


Coleman-Weinberg potential at zero temperature [44]. 
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where 7 runs over all particles in the NMSSM, with each having degrees of freedom ni, 
field-dependent mass m;(y;) and spin s;. In this work, we adopt the Landau gauge and the 
DR scheme with renormalization scale Q to calculate Vow. In order to maintain tree-level 
relations involving VEVs, counter terms Vor should be introduced. Vr denotes the thermal 


correction at the finite temperature, 
-5 2 in f dx x? nU Fe VE t/t), (8) 


which can be improved by daisy resummation contributions Vagisy. In practice, this term 


can be achieved by inserting thermal mass contributions in the field-dependent mass. 


To avoid the baryon asymmetry generated at the EWPT being washed out, the phase 
transition must be strongly first order, which can be quantitatively characterized in the 
perturbative calculation as 


ve(Te)/Te 2 0.9. (9) 


Here T, and ve are the critical temperature and order parameter of the phase transition. 

We implement the effective potential Eq. [6]into the public package CosmoTransition [46] 
to evaluate the phase transition numerically and numerically perform a parameter space 
scan. As pointed out in previous works [45], there are mainly three phase transition 
patterns in the NMSSM. When the temperature decreases, the Universe can directly transit 
from the symmetry phase into the electroweak breaking vacuum Qpw (Type-III), or undergo 
a intermediate phase in the singlet subspace (Type-I) or the H, subspace (Type-II). It has 
been found that a strongly phase transition can be achieved in Type-I and Type-III transition 
without light squark contributions to thermal loops [42]. 

The shape of the effective potential at zero temperature, more exactly, the energy gap 
between the symmetry phase and the broken phase (AV = Vgym — Vew) |r=o encodes the 
information on EWPT in the NMSSM. A smaller AV may lead to a lower T, and thus a 


large transition strength ve/Te. This correlation can be understood by 


T. is G Ir-r.) AV (10) 


Following this insight, we can perform a semi-analytical analysis of the transition strength 
in terms of the energy gap [42]. 
First of all, the potential energy of the electroweak vacuum can be divided into three 


parts at tree level 


Vew = Viy + Vidw + Vile, (11) 
where the individual contributions are given by 
y2 dM f vem? 
Van = = 6 (cos? 26 + J sin? 28) a=] P, (12) 
1 
wre — 3h AKUs —K'v,, (13) 
A 
Ve = — (1 = 7 sin 28 — = sin 28) pv" = —Cap?v" i (14) 


Fa 


Ld 


Explicitly, the first part is the contribution from the doublet and is almost fixed by the SM- 
like Higgs mass. The second part is the pure contribution from the singlet. The third part 
results from the doublet-singlet mixing, which can be described by an auxiliary parameter 
Ca, 

A) 


Ca =1— (Gè + 5) sin 26. (15) 


In the case of type-I phase transition, the universe transits into an intermediate phase in 
the singlet subspace before the EWPT. The absolute minimum of the singlet subspace at 


zero temperature locates at the origin or 


—A, 
Us = Te (1+ vi- 82, ). (16) 
where x, is given by 
1 1 
Tk = 8 B gil + Ro z Car?v? /A2, (17) 


where R, = 4k/AA,. The Tpyp-I transition occurs only if the us is the absolute minimum 


with xz, < 1/9. The relevant energy gap can be easily obtained as 


AV = Vs — Viy — Viw — Vee 


22 
we h — Cav (u? — o?) + 6? (0? — u?)? 

1 
+ 3 tAn [2u (us = Us) + vsv? = uz) ’ (18) 


where Vg is the tree-level potential energy of us. It is obvious to see that a substantial 
deviation of u, from v, is crucial to decrease the energy gap away from the doublet limit 
set by the Higgs mass. The numerical results showed that the Type-I phase transition 
correspond two cases, in which either Rẹ < 10 and —1 < Rẹ < 0 [42]. For the large |R], 
the deviation of u, from v, would be small and is disfavored by the SFOEWPT. 

Note that at tree level the vacua energy gap AV may be negative. Fortunately, the 
Coleman-Weinberg one loop correction would lift Vs and reduce Vew, thus guarantee that 
the EW vacuum Qpgw is lower than Qs. However, if the metastable vacuum Qg is well below 
Qew, the finite temperature potential could not smooth out a large negative gap when the 


Universe cools down. In this case, the phase transition is not valid. 


Then we consider the Type-III phase transition arising in the case that the origin is the 
(19) 


absolute minimum in the singlet subspace or the origin is metastable. In this scenario, the 


+1). 


tree-level energy gap is simply given by —Vew as 
2 
um), 2,2 2,4 
C 
+ Cau v +K ETA 


AV ~ 
4 
We can see that the singlet part will dominate the energy gap for a large u. A small energy 


gap often requires a large negative A, and —4/3 < R, < 0. For a moderate ju, the mixing 


part becomes important and a negative C4 is useful to decrease the gap. 


Ill. NMSSM INTERPRETATIONS OF GALACTIC CENTER EXCESS 
The Galactic center excess( GCE) prefers a ~30-40 GeV DM particle annihilating directly 


into bb with a cross-section about (av) ~ 2 x 10776 cm3/s. In this work, we would like to 


explain this excess by the annihilation of neutralino pair to bb via an s-channel light CP-odd 


pseudo-scalar in the NMSSM. 
The neutralino mass matrix in the NMSSM is written as 
givu 0 


a E 
M, Sa ne 9 
M= 0 =u Àv (20) 
0 —\va 


The lightest mass eigenstate of the neutralino is the DM candidate, which is defined as: 
(21) 


XS NyB + N W + Nia fa + Nui, + N58. 


The Majorana DM annihilation cross section to bb with a relative velocity v, for the 
(22) 


3 


interaction $a(YayyXYX + Yarwby’b) is given by 
P Yx Yabo 
32rm? (6 + v2/4) + 9?’ 


OU, 


2 


where y = Mala /4ms, 


is a parameter defined by the pseudoscalar decay width Ty, 6 is a 


degeneracy parameter defined as 


m2 


= — ae 2 
5=1 ima (23) 


If 6 is not very small, the DM annihilation occur off-resonance and is almost velocity inde- 
pendent. As discussed in Ref. |8], the GCE and correct thermal DM relic density can be 
simultaneously accommodated as long as the product eae is adjusted to an appropriate 
value. In this case, the main components of the DM would be Bino and Higgsino in the limit 
k/A > 1. Since the singlet scalar is heavy and is almost decoupled, the SM-like Higgs is 
similar to that in the MSSM. The singlet component of the lightest pseudoscalar can sup- 
press its production cross section, hence the stringent limits on the the pseudoscalar mass 
from the LHC H/A — r*7~ searches can be avoided. However, large «/ would enhance 
the -Vw in Eq. (13), which is disfavored by the SFOEWPT. Furthermore, a very small A 
is difficult to give a correct SM-like Higgs mass ~ 125 GeV at tree level. Therefore, in this 
work we do not consider such parameter space. 

In the limit K/A < 1, the dominant component of the DM may be Singlino. In this case, 
the correlation between «/A and the DM mass m, in the Zs NMSSM is 


Km |, Av? (m, — usin 28) My 


24 
à 2p my(m? — u?) 2u’ a 


where the second equality corresponds to small m,/u. The DM would also have non- 
negligible Higgsino components with a moderate u value, which can be described by for 


small my /u 


N 

a YY gg BG an) (25) 
Nis H u 

Nia AU i My 

— ~ —— sin f(1— . 26 
Nis H l T oe 


For the CP-odd Higgs, the absence of tachyon states favors a negative KA, or a small 
positive KA,. The constraints from the LHC Higgs researches also compress the allowed 


parameter space. If the heavier CP-odd Higgs is very heavy which is consistent with the 


10 


Higgs observations, the lightest CP-odd Higgs mass can be written as 
1 AA 
m? ~ zA sin 28(7 + 4K) — BTR Ay, (27) 
u 


If Av sin 23/2, is not very large, the main component of the lightest pseudoscalar is singlet. 
In this case, the coupling gabb is suppressed by a small active pseudoscalar fraction. At the 
same time, the coupling gay, is determined by small x and small Higgsino fraction. This 
means present DM annihilation cross section should be enhanced by the resonance effect 
with small 6. From Eq. we can see that the GCE can be easily explained by adjusting 
the combination gayx9abb/ (0° + 77) in the zero temperature limit. In the early Universe, the 
DM annihilation occurring close to the resonance is sensitive to the temperature [£7 [48]. 
If 6 < 0 and |6| is not very small, the øv, tends to have a larger value with a larger v, as 
the temperature increases in the early Universe. In this case, the DM relic density would be 
suppressed. Some non-thermal DM production mechanisms are needed to obtain the correct 
DM relic density. We do not consider such mechanisms and focuss on the thermal freeze-out 
mechanism here. For the 6 > 0 case, from Eq. we can see that the ov, always decreases 
as the temperature increases. This means the process yy — a — bb can not sufficiently 
reduce the DM abundance in the freeze-out epoch. Some other DM annihilation channels 


are needed to generate the correct DM relic density. 


IV. NUMERICAL RESULTS 


In this section, we perform a numerical analysis to interpret the GCE and SFOEWPT 
in the NMSSM. We use the packages NMSSMTools 4.4.0 [49,50], micrOMEGAs 3.6.9.2 
and CosmoTransition [46] to scan a six-dimensional parameter space (A, K, Ay, Ax; 4, tan 8) 
with various experimental limits. Since we focus on the EWPT properties affected by the 
singlet sector, we assume that all the sfermions are heavy and fix soft breaking parameters 
as Mr = Mg = 200 GeV and Mg = M, = Ma = At = Ap = A; = 2000 GeV. A benchmark 
point, which can simultaneously explain the GCE and SFOEWPT, is given in Tab 


In order to understand the properties of the parameter distribution favored by the GCE 
and SFOEWPT, we plot Fig.|1} where only two parameters « and A, vary in the parameter 
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À A (GeV) & A, (GeV) tang u (GeV) 


0.5 840.0 |0.029 -99.18 3.05 235.0 


mgo (GeV)| Mma (GeV) QR? |(ov)|y-0(cem?/s)| ogy(cm?) PTS 


34.96 69.54 |0.095) 1.72 x 107-26 |9.03 x 1071?! 1.06 
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4 


0.022 30 0.032 0.034 0.022 0.024 0.026 0.028 0.030 0.032 0.03 
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FIG. 1: Results of a numerical scan with fixed: A = 0.5, A, = 840 GeV, u = 235 GeV, tan 8 = 3.05, 


Mı = M2/2 = 100 GeV, M3 = 800 GeV, Mr = Mp = 200 GeV, Mg = My = Ma = At 

A, = A; = 2000 GeV . The green region denote the parameter space where the DM annihilation 
ov, ~ (0.5 — 4) x 10776 cm? s~! can interpret the GCE. In the left panel, The blue region can 
explain the correct DM relic density 0.091 < Qh? < 0.138 via the Z mediated annihilation 
in the early Universe. Also shown are contours of ve/Te (the magenta dotted lines) which is usually 
required to be larger than 0.9 by the SFOEWPT. The red and gray regions are excluded by the 
direct detection LUX and the Higss mass 124 GeV< mp < 128 GeV. In the right panel, we 
show contours of the DM mass my (red dotted lines) and mass difference 2m, — Ma (blue dotted 


lines). 


space near the benchmark point and the other parameters are fixed. In the right panel of Fig 
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FIG. 2: Results of a numerical scan with varying A and u, k = 0.029, A, = —99.18 GeV and other 


parameters being fixed as Fig.|1| Other discerption are same as Fig. 


red and blue lines denote DM mass m, and mass difference of 2m, — Mma. The behavior 
of these lines can be easily understood from Eq. and Eq. In Fig [1] we show the 
region consistent with the GCE with an annihilation cross section gv,|,,-39 in the range of 
(0.5 ~ 4) x 107% cm? s~!. As expected from Eq. there are two distinct strips near the 
contour 2m, — Ma = 0 which correspond to the cases ô > 0 and ô < 0 respectively. 

If there is no other DM annihilation channel, it is difficult to explain the GCE and correct 
DM relic density simultaneously. In the lower green band for 6 < 0, the DM annihilation 
occurring closer to the resonance in the early Universe would induce a very small relic density, 
while in the upper green band for 6 > 0, the DM particles would be overproduced due to 
the suppression of the DM annihilation cross section in the early Universe. In the left panel 
of Fig !1| the blue region denotes the parameter space where a correct DM relic density can 
be obtained. In this region, the DM annihilation via the exchange of the Z boson plays an 
important role in the early Universe, and sufficiently reduce the overproduction of DM for the 
ô > 0 case. This annihilation channel is a p-wave process and can be negligible in the Galaxy 
with small DM relative velocity. The green and blue bands can overlap at (K, A,) ~ (0.029, 
-99.2 GeV). Since the correct DM relic density requires a non negligible DM-Z interaction, 
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we also check the partial decay width of the Z to a pair of DM which is controlled by the 
Higgsino component of DM. For our benchmark point, this width is 0.59 MeV and is allowed 
by the upper limit ~ 2 MeV from the Z decay measurements [54]. 

Then we show the parameter distribution of (A, u) favored by DM results in Fig |2| where 
the other parameters are taken in Tab JI| Similar to Fig !1| there are two bands for the GCE 
correspond to different ô values. Eqs. and [26]imply that the Higgsino faction of the DM 
and thus the DM coupling to the Z boson are controlled by A/u. By tuning the value of 
A/p, the correct DM relic density can be obtained via the Z channel annihilation in the blue 
region. The upper green band and blue region overlap at (A, u) ~ (0.50, 235 GeV). 

In Fig. and |2| we also show the constraint from the direct detection. The LUX 
collaboration has set stringent constraints on the DM with mass ~ O(10)GeV [55]. The 
spin-independent signature for our benchmark point is dominantly induced by the DM- 
quark scattering process via the exchange of CP-even Higgs. The scattering cross section 
is strongly affected by Higgsino components in the DM. This explains why the parameter 
space with large À and small u would be excluded by the LUX results. 

We depict the contours of v./T, in Figs |1} and It is can be found that there does 
exist a parameter space satisfying the DM results and SFOEWPT. Eq. implicates the 
relation between the phase transition strength and the energy gap, deeper understanding of 
the SFOEWPT calls for a detail analysis of the energy gap. The energy gap at tree level 
AV iree and full numerical result with the Coleman-Weinberg one loop correction AVpu are 
shown in Fig. |3| We find that the variation tendency of phase transition strength ve/Te with 
respect to A(x) and ju(A,,) is consistent with those of the energy gaps AV;ree and AV fu. 

In our benchmark scenario, Type III SFOEWPT is preferred since the mixing term with 
a negative C4 tends to sufficiently decrease the energy gap. The parameter dependence of 
the AV;ree can be easily understood by Eq. [19] which is rewritten as 


22 
v mi, 


AV ~ 


1. k. 1 
— 5sin 2BAyuv + (1 — Ë sin20)pu? + TACT af mE) (28) 


where the first term, middle two terms, and last two terms represent the contributions 
from the doublet, mixing, and singlet parts, respectively. For our benchmark model with 


k ~ O(0.01) and A, ~ —O(100) GeV, a negative contribution of «A,v2?/3 dominates in the 
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singlet part. Hence, AV;ree decreases as « and |A,| increase in the region of k < —A,/6u; ~ 
0.035, and then the phase transition is enhanced as shown in Fig. ]and Fig. [I] From Eq. 
we also find that AV;ree is suppressed by the small À in the mixing and singlet parts. The 
dependence of AV;iree on u is more complicated. In our benchmark scenario with a small x, 
the contribution from the mixing part would be more important than that from the singlet 
part. In this case, it can be found that AV,,.. becomes more negative as u decreases when 
u > Aysin26/4 ~ 110 GeV. 

Note that although the AV at tree level is negative in Fig. |3| the loop corrections from 
the Coleman-Weinberg potential lifts up the AV;,-. and helps the phase transition to occur 


in the early Universe. The detailed discussions of the loop corrections to AV can be found 


in Ref. [42]. 
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FIG. 3: The energy gaps at tree-level AV;ree ( black contours) and the full numerical results 
with loop corrections AV+,,, (magenta contours). Left panel: plots of AV in the k-Aẹ« plane with 


parameters set as those in Fig.|1} Right panel: plots of AV in the A-y plane with parameters set 


as those in Fig. 
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V. CONCLUSIONS AND DISCUSSIONS 


In this work, we study the possibility that the GCE and SFOEWPT can be simultaneously 
explained in the NMSSM model. The GCE can be interpreted by the annihilation process 
xx — bb via the exchange of a light pseudoscalar. For the Singlino/Higssino DM, the 
annihilation should occur near the resonance with a very small mass difference 2m, —m, > 0. 
In this case, the correct DM relic density is obtained by the DM annihilation via the exchange 
of a Z boson in the early Universe, meanwhile the present DM annihilation cross section 
accounting for the GCE can also be easily achieved. 

Besides the mostly singlet pseudoscalar ~ O(10) GeV required by the GCE, there is 
also an accompanying light CP-even Higgs in the scalar sector. Such singlet like scalar 
could decrease the vacua energy gap AV with a negative A, and a moderate u, which 
would lead to SFOEWPT required by successful EWBG. We find that there does exit the 
parameter space in which the GCE and the correct DM relic density can be interpreted by 
the Singlino/Higssino DM with SFOEWPT being realized at the same time. 

The discussions in this work can be extended to the DM models containing new singlet 
states. In such models, the Higgs spectra would be affected by the requirements from the 
EWPT and DM phenomenology. Discovering such new scalar sector at O(10)GeV is an 


excellent motivation for the future high energy collider experiments. 
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